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An elementary treatment of the Dirac equation in the presence of a three dimen- 
' sional spherically symmetric delta potential is presented. We show how to calculate 



o 



the cross section using the relativistic wave expansion method for a one delta po- 



\ tential and two concentric delta potentials. We compare our results with the cross 

section calculated in the Born approximation. 



> . I. INTRODUCTION 

o 

cn 

^ ' The non relativistic scattering theory, and in particular the partial wave expansion 

^ , method, is a very well known issue. An extensive literature ^j] exists in this subject. In 
O 

discussed. This has to do, probably, with the fact that relativistic quantum mechanics is, in 



the relativistic case, however, the scattering problem with a potential, has almost not been 



^ \ fact, a theory of many particles and, therefore, quantum field theory is a more appropriate 



language to discuss these kind of problems. 

We think, however, that the scattering problem in the frame of the Dirac equation deserves 
^ ' some attention. This fact has motivated us to make use of a simple mathematical method, 
valid for central potentials, the so called relativistic partial wave expansion method provides 
us with valuable information for relativistic particles near the low energy limit, giving us 
also values for the phase shifts and cross section of the scattered wave for different channels 
of the angular momentum. This method will be used to complement the work done in Q], a 
discussion for bound states in the presence of a delta potential, which shows how to fix the 
boundary conditions for the wave function when crossing the V{r) = a6{r — tq) potential. 
We will like to remark that there are general treatments for boundary conditions for this 
kind of singular potentials, in the frame of the construction of self adjoint extensions for the 
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Dirac Hamiltonian P|. 

In this work we use the relativistic partial wave expansion method with the boundary 
conditions mentioned above, making a simple analysis for the scattering states. Moreover, 
we will also introduce a double delta potential of the form V{r) = ±ai6{r — ri) ±a2S{r — r2). 

The idea is to discuss the behavior of the differential and total cross section when we 
vary the energy of the incident particles, analyzing the occurrence of certain resonances in 
different channels. We also find a useful relationship between the phase shifts and the total 
cross section. 

This work is complemented with the use of the Born approximation, that gives us a more 
general perspective of this problem. 



II. RELATIVISTIC SCATTERING 

As it is extremely well known, asymptotically, we expect that the wave function in the 
non relativistic case will behave as 

<i,{r)r^^exp{tKz) + f{e,<P)^^^^^. (1) 

where K is the momentum. The differential cross section is then 

Ho,<P) = ^ = \m<i>)\'- (2) 

For scattering at high energies of 1/2-spin particles, we must consider the Dirac equation 

{ca ■ p + (3mc^ + V{r))ilj{r) = Hilj{r), (3) 

where a and (3 are the usual 4x4 Dirac matrices, and is the four- component Dirac spinor. 
We know that for a free particle, the one dimensional solution of the above equation is 



where e = ±Ep, a = ±1, Xi 







2e \ s^P—y 
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When working with a spherical symmetric potential it is better to transform the Dirac 
equation into spherical coordinates 0,0], where the complete set of commuting operators is 
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given by if, J^, J3 and k, where J is the total angular momentum and k is defined by 



p{-S ■ L + hl). 

f V 

The eigenvalues of the operator n are given by 



k^ = = ±{j + 



(5) 



(6) 



If we parametrize the four-component spinor by separating the radial and angular depen- 
dence according to the Ansatz 



^ 9{r) XU9. 0) ^ 



where 



(7) 



(8) 



being C{l^j; fi — m;m) are the appropriate Glebsch-Gordan coefficients. It is easy to see that 
the k operator satisfies the following relationship with the angular momentum-spin function 



(9) 



For the two different coupling, j = I + 1/2 and j = I — 1/2, the eigenvalues are given by 

(10) 



I j = Z - 1/2 K > 

-/ - 1 j = 1 + 1/2 K<0 



The dependence of k on the eigenvalues of the angular momentum will be extremely 
important for our calculations. Everything will be function of k. However, in the future 
we will replace k with the respective 1^ remembering the existence of the LS coupling. We 
notice that: 



and 



K K > 

-K-1 K <0 



K-1 K> 
—K K < 



(11) 



(12) 
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The Dirac equation in spherical coordinates has the form 

(^-ic&r + ^ - 7'^) + ^"^^^ + '^(^)) ^(^) = ^^(^)- (13) 
The solution of this equation for a free particle with the Ansatz given in ((7j) is: 

G(r) = r(aix(r) + hm^ir)), (14) 

= |^ ^^^^^2 ^("i-^'^->c(^) + bm^.i^)). (15) 

where G'(r) = rg{r), F{r) = rf{r) and ji, rii are the regular and irregular Bessel functions 
respectively. 



III. PARTIAL WAVE EXPANSION METHOD 

Here we will present the relativistic wave expansion method. The reader can see the 
details in -S]. This method shows us how to calculate the differential and total cross section 
just by knowing the phase shift of the scattered wave. 

Let us consider the asymptotic behavior of a 1/2-spin scattered wave 

^,^ = a,e^v{^Kz)+h{e,ct^Y-^^ ^ = 1,...,4. (16) 

To be able to find the effects of the scattering in the spin of the particle, we must 
first notice that not all the components of the spinor that represents the free particle are 
independent, which, in turns, allows us to see that not all the bi{6,(j)) components are 
independent. For example if we have a spin parallel to the direction of incidence (spin up), 
we have: 

^i.^ooexp(^irz) + /+(0,0)^^^^^, (17) 

^^J,^^^g+(9,^)^^^^^. (18) 
Now, if the spin is anti-parallel to the direction of incidence (spin down) we have 

^.^..-(^,0)^^^, (19) 

^^^^^exp(iKz) + r(e,(j))^^^^^^^. (20) 

r 

The functions /^(6', 0) and g^{0,(l)) are called scattering amplitudes (similar to the non 
relativistic case). With these amplitudes we can calculate the differential cross section 



5 



for a non polarized beam, making the change = f = f, 9^ = 5'exp(i^) and g = 
—gexp{—i(f)) getting 

m<P) = % = \f\' + \9\^ (21) 

where the new scattering amphtudes are given by 

oo 

f{e) = Y.AiPi{cos9), (22) 

1=0 

OD 

gi9) = Y.BiPlicos9), (23) 
1=1 

and 

= + l)[e^P (2^^/") - 1] + ^[e^P (2^^r) - 1]}, (24) 

= ^[^^P ^^zr^r) - exp {2ir]+)]. (25) 
We can see that the differential cross section has a direct dependence on the phase shifts 
of the wave function. The notation rif" gives us the dependence of the wave function on the 
angular momentum and on the two different coupling (j = / ± 1/2). It is easy to calculate 
the total cross section, just by integrating the differential cross section, obtaining 

'^totai = 4vr ^ . (26) 

IV. 1-DELTA POTENTIAL 

As we mentioned before, we will apply the wave expansion method to a spherical potential 
given by 

V{r) = ±a6{r -ro). (27) 

To be able to find the phase shifts, in order to obtain the differential and total cross 
section, we .us. use the boundary condition. desCbed in Q: 

Fl + Gl = F^+G^, (28) 

and 

F\ ^ (F./G.) + a 

G+ l-a(F_/G_)' ^ ^ 

where F+ y are the solutions of the Dirac equation outside the potential (r > ro), and 

-F_ y G^ are the solutions inside the potential (r < tq). a is an dimensionless constant that 

involves the constant a, a = tan [a /Tic). 
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The first boundary condition tells us that the absolute value of a function with real part 
F and imaginary part G is constant, which is nothing but the continuity of the probability 
density. The second boundary condition will be essential to calculate the phase shifts. 

First we must find the solution of the Dirac equation in all regions of space. Since we are 
dealing with a delta type potential, it will only influence the wave function in the support 
r = tq, being, therefore, the solution is the same as in the free particle case. Separating the 
space in two regions: 

I. - For r < Tq : 

Gi{r)=rarjakr), (30) 
\k\ E + mc^ 

The functions ni^{kr) and ni_^{kr) do not appear in the above equation, because the 
wave function must remain finite at the origin. 

II. - For r > ro 

G2{r) = r{ji^{kr) cos{r]f) - ni^{kr)sen{r]f)), (32) 

F2{r) = YITT-, 2'^ji^-'^kr) cos{r]^) - nu^sen{rif)). (33) 

Here we have written a2 and 62 as the phase shifts of the scattered wave. These phase 
shifts can be found using the boundary conditions. Equation ()28|1 does not give us much 
valuable information. It only allows us to find the constant ai in the wave functions. We 
will concentrate on the second boundary condition (j^ . Here the constant ai disappear, 
leaving us 'qf' as the only variable. A simple algebra gives us: 

tan(?7r'') = — / , "^'".^ 775 : — r, (34) 

where 

^=^_^. (35) 

\k\ E + mc^ 

The spherical Bessel functions can be obtained through the following recurrence relations 



a 
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FIG. 1: Behaviour of the first three phase shifts for different values of the energy in the one delta 
potential case (in units of mc^), r/o (solid curve), r/j~ (dotted curve), -qf (dashed curve). The super 
± index corresponds to the different angular momentum-spin couplings, the sub index indicates 
the angular momentum channel. Here we use a = —Ihc and r = fi/mc 



The undulatory behavior for the phase shifts in Figure 1 shows the typical dependence 
on the Bessel functions. Notice the occurrence of discontinuities in each angular momentum 
channel, for certain value of energy. This discontinuities, which we will identify as resonances, 
appear when the -qf get the value n7r/2, with n a natural number. In our work we will only 
discuss the case where n = 1,-1. However, we will see that this resonances in the phase 
shifts, will not affect the cross section. 

Now, using the results of the phase shifts, we can find an expression for the differential 
and total cross section using the equations (pT|l and for different values of the angular 
momentum. 

Figure 2 shows the total cross section for different values of the energy. Here Z = 1 is the 
sum of / = 0, 1 and / = 2 is the sum of / = 0, 1, 2 (see equation ^^). 

We first notice that those resonances in the phase shift analysis do not seem seem to have 
any consequences in the total cross section, in spite on their direct influence in equation 
(I2(i|l . This has to do with the fact that the resonances in the phase shifts appear at different 
values of energies. The cross section vanishes smoothly for high energies. For low energies 
(near mc^) the dominant channel of the angular momenta is / = (as expected). It is also 
important to mention that if we increase the delta radius or the coupling constant a from the 
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1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 



Energy/mc 

FIG. 2: Behaviour of the total cross section for different values of the energy in the one delta 
potential case (in units of mc^), axil = 0) (solid curve), axil = 1) (dashed curve), (Tt(^ = 2) 
(dotted curve). Here we use a = —IHc and r = H/mc 



potential, the total cross section decreases and vanishes faster when we increase the energy. 



V. 2-DELTA POTENTIAL 



Let us consider a two concentric delta potential: 

V{r) = ±ai5{r -ri)±a2S{r -r2). (38) 

We will carry out the same procedure as before. We separate the space into three regions, 
finding the corresponding solution of the Dirac equation. Our potential will only affect us 
when r — ri,r2 and, hence, the solutions are the same as for a free particle: 

1) For r < ri 

Gr{r)^ra,jiAkr), (39) 
\k\ ti + m& 



2) For ri < r < r2 



G^ir) = r{a2ji^{kr) + b2ni^{kr)), (41) 
= 0^^-^(«2ji_.(^r) + b2ni_^{kr)). (42) 



9 



3) For r > r2 

G^{r) = r(ji^(kr) cosr]f — ni^{kr) smrff), (43) 
-^3(^) = TlT^~r~~^Ui'^(^^) cosrif" - ni_^{kr) smrif). (44) 

t^j Hi ~r TflO 

Using the boundary conditions for each deha and having in mind that a\ = tan (ai/hc) 
and ^2 = tan {a2/hc), we can find an expression for the phase shifts: 

«2 (A^jI^ + jI + A^ni_ji_^ + mji) + A{kri) (A{niji_^ - ni_jij) 

tanr^f = — ^ ^ ^ (, (45) 

M^i-Ji. - niji^j + a2 [A'^ni^Ji^^ + riiji^ + A{kri){A^nf_^ + nfjj 

where 

K khc 

A 



\k,\ E + mc^ ' 

and 

ai (A^if + i? , 

^(^^i) = ■ .\ % (46) 

In equation ()45|1 all Bessel functions depend only on kr2- The dependence on kri is 
concentrated on the function A{kri) as shown on (j46|) . Now, if we choose Oi = or ri = r2 
the problem, and hence the solution, reduces to a one delta potential. 

Knowing the phase shifts for a two-delta potential, we can obtain an expression for the 
total cross section. Figure 3 shows the cross section as a function of the energy. 

Again we can see that the resonances of the phase shifts do not have any notorious 
infiuence on the total cross section, in the same way as happened for a one delta potential. 
It decays smoothly to zero for high energies. For small energies, near mc^, the dominant 
angular momentum is the I = channel. Once again, we must remember that I = 1 stands 
for the sum of / = 0, 1, and / = 2 stands for the sum of / = 0, 1, 2. 

VI. BORN APPROXIMATION 

Let us remember our potential 

V{r) = ±ai6{r — ri) ± a26{r — ri). (47) 



Following the same steps shown in Itzykson-Zuber |^ (equation 2-126), we get for the 
scattering amplitude: 

Sf^ = — ^^(2vr)5(E; - E,) f d\exp{-zq- f)V {r)u" {p fh%^ (p,) , (48) 
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1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 



Energy/mc 

FIG. 3: Behaviour of the total cross section for different values of the energy in the two delta 
potential case (in units of mc^), aril = 0) (solid curve), axU = 1) (dashed curve). Here we use 
ai = Ihc, 02 = Ihc, r\ = 2h/mc and r2 = 3h/mc 



where q — Pf — Pi- Then r • r— \q\\r\ cos6 — \q\\r\^, with = cos^. 

Integrating in spherical coordinates: J^d/j, Jq^ d4> J^r'^dr, defining F = u"{pf)'~^^u^{pi), 
we obtain 



mi 



S{Ef — Ei)[±airi sin qri ± a2r2 sin qr2\T. 



V{EiEf)-2 q 

The transition probability from i to f per time and particle is: 



(49) 



dt 



I 



^^'^{±a,n sin qn ± a2r2 sin gr2)r 



x2nS{Ef-E,)V^. 



A straightforward calculation gives us 



dt 



m 



V{E,Ef) 



{±a,nsmqn ± a2r2sii^|^|,^^^^ _ 



k1 



(50) 



(51) 



To be able to find the differential cross section, we must divide the last equation by the 



incident flux 



- — (±airi singTi ± a2r2 sin gr2)^|r|^(5(£^/ - Ei)p'jdpfdfl. 

\Pi\-^f\q\ 

Using the fact that \pi\ = \pf\ = Pf and Pfdpf = EfdEf we obtain 



da 

dn 



m 



(±airi sin qri ± a2r2 sin gr2)^|r| 



(52) 



(53) 
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Since we are interested in the cross section for a non polarized beam, we must sum over 
all a and (5 that appear on the function F, this means 

di(J 771/^ 1 

— = Y^{±airismqri ± a2r2smqr2fY.-^\u'^{pf)-f^u'^{pi)\'^. (54) 

We will concentrate directly on the sum over a and p. Rewriting the sum as traces, we 
get 

E 2 E l« (P/)7 U^{P^)\ = -^tr (^7 2m ^ 2m J ' ^^^^ 
and using some identities for the traces of the 7 matrices, we have 

tr{^\lW{l'pf) = -Pi-Pf + E,Ef), (56) 

tr(707°) = 4. (57) 
We will also need some kinematics relations 

Ei^Ef + E, (58) 

n 

Pi-pf^E^- p^ cos ^ = m^ + 2£;2/?2 sin^ -, (59) 
where (5 = v/c= \p\/ E is the incoming velocity. In this way we obtain 

,ao) 

The differential cross section with \q\'^ — 4|p|^sin^ | = 4/3£'|p| sin^ | is: 

±a2r2sin(2y^^sin|r2)]2. (61) 

We can easily see that the differential cross section is the same for attractive or repulsive 
potentials. One important thing to notice, when working with the partial wave method, we 
focus on different channels of the angular momentum. But, however, we take the sum of all 
the angular momenta in the Born approximation. 

The total cross section is given by 

(^totai = j d9sm9d(f)^, (62) 
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where da/dfl is defined in dHD), We obtain the following equation: 

2ttE 6 9 

(ytotai = — r-T / d9 sin 6'[air^ sin^ (ai sin -) + ajr^ sin^ (aa sin -) 
4:p\p\ Jo 2 2 



. . . . . . / I - r sin^f \ 

+2aia2rir2 sm [ai sm -) sm (^2 sm -)J — ^ . (63) 

2 2 V sm 2 / 



where ai = 2^ [3E\fyri y 0:2 = 2y (3E\pr2. We must notice that the factor ± is not taken 
into account. As before we realized that there will be no influence on the results if both 
delta potentials are repulsive or attractive. The case where the sign of ai is different from 
the sign of a2 will also not give new information. Integrating directly (jUHj) we obtain: 

ttE 

(Ttotai = — -{-2«(2ai)(airi)2-2«(2a2)(a2r2)' + 21n(2ai)(airi)2 
2(3\p\ 

+2 In (2a2)(a2r2)' - f3\{ainy + (a2r2)') 
„2 /sin(2ai), , .n sin (2a;2) , .o\ 

V «1 "2 / 



^9 /sin^ (ai) , ,2 sin^ (q;2) 
^ ^ am ' + 

V «1 «2 



+A{airi){a2r2){Ci{ai - ^2) - Ci{ai + «2) 
^2rSin (ai - 0:2) _^ cos (ai - ^2) 1 

«! — ^2) («!— ^2)^ («! — «2)^ 

sin («! + ^2) cos (ai + ^2) , 1 



(ai + a2) {oLi + a2Y (ai + 0:2)^ 
-C(2(airi)2 + 2(a2r2)2), (64) 

where Ci is the well known cosine integral, given by the following equation: 

«(x) = C + In (x) + E(-l)'^^' (65) 

and C is the Euler constant. 

Figure 4 shows the dependence of the total cross section while we vary the energy of the 
particles. We notice that the total cross section has the same behavior found in the partial 
wave method. On the other hand if we take a Ar = r2 — ri = cte, we can see that the 
behavior is again similar to the one we found in the partial wave method. 

We can say that the behavior of the cross section given by the sum of all the angular 
momenta is similar to the behavior of the cross section for a specific angular momentum 
channel. In spite of several efforts, we were unable to find geometric conditions related to 
the separation of the two delta potentials, and the possible occurrence of resonances as peaks 
in the total cross section. 
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FIG. 4: Behaviour of the total cross section for different values of the energy in the two delta 
potential case (in units of m(?) with the Born Approximation. Here we use oi = —Ihc, 02 = —Ihc, 
ri = 2h/mc and r2 = 3h/mc 



VII. CONCLUSIONS 

The general idea of this work was to use the boundary conditions obtained in Q], when 
working with delta type singular potentials, in particular for a spherical potential like V{r) = 
±a6{r — tq) in the scattering region. We studied the behavior of the phase shifts, the 
differential and total cross section. 

We found a direct relation between the phase shifts of the scattered wave and the cross 
section. With this we were able to calculate values of the cross section for different angular 
momentum channels, giving us valuable information, in particular, the dependence of the 
cross section on the / = channel in the low energy limit. 

It is important to notice that the behavior of the cross sections for each channel of the 
angular momenta is very similar for repulsive and attractive potentials, in the one and two 
delta cases. We confirmed our results through the calculation of the cross section within 
the Born approximation. This method does not require the knowledge of the boundary 
conditions, but it does not give us information about the angular momenta of the particles. 

It could be interesting to expand this work to an anion-potential, i.e. a series of concentric 
delta-potentials, trying to find recurrence relations for the occurrence of resonances in the 
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different angular momentum channels. 
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